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*S3 RELATED APPLICATIONS 

This application is a continuation-in-part 
yS application of Application Serial No. 09/124,339 filed 

ff; July 29, 1998, entitled, "Method and System for Modeling, 

is i 5 Predicting and Optimizing Chemical Mechanical Polishing 

s Pad Wear and Extending Pad Life". 

5§ TECHNICAL FIELD OF THE INVENTION 

4; The present invention relates generally to systems 

Q 10 for the chemical mechanical polishing (CMP) of workpieces 

such as semiconductor wafers. More particularly, the 
present invention relates to a pseudo-physical model of 
both the contact forces between the polishing pad and the 
wafer and the erosion of the wafer during the polishing 
15 process. Specifically, the pseudo-physical model of the 
present invention simulates the result of a CMP process 
predict eroding topography and thereby facilitate the 
optimization of the process parameters associated with 
the operation of an actual CMP system. 
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BACKGROUND OF THE INVENTION 

Chemical mechanical polishing (CMP) of semiconductor 
wafers has become the preferred method for planarizing 
dielectric or metallic layers at various stages of 
integrated circuit fabrication. During the CMP process, 
a workpiece surface, such as the surface of a 
semiconductor wafer, is held against a rotating platen. 
The rotating platen is covered with one or more- 
slurry- soaked polishing pads. The combination of the 
chemical interaction of the slurry with the semiconductor 
wafer and the friction between the polishing pads and the 
surface of the semiconductor wafer removes material from 
the wafer. Because of the small dimensions of silicon 
wafers and the precise surface finishes required and the 
costs associated with changing polishing pads, there is a 
need to gain a better understanding of the CMP process by 
mathematically predicting the results of a chemical 
mechanical polishing process. 

Although the bulk material removal rate of 
interlevel dielectric materials by a CMP system is 
heavily chemistry dependent, there is no chemical 
preference between the high and low areas on the surface 
of the wafer being polished. Therefore, when using a CMP 
process on dielectric materials, the planarization of the 
polished surface of the wafer is due solely to the 
mechanical action of the polishing pad against the 
surface of the wafer. And it is this mechanical action 
which is predictable by the use of pseudo-physical 
modeling techniques . 

CMP systems are often used to process the silicon 
dioxide dielectric material (commonly referred to as 
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"oxide") layers from semiconductor wafers. The objective 
of the use of such CMP systems is to uniformly remove the 
dielectric material or oxide across the semiconductor 
wafer being polished such that the small (sub-micron to 
5 millimeter) features or surface irregularities that 
populate the wafer surface are worn away and thereby 
eliminated. It is important, at the same time, to assure 
that, the overall global characteristics of the wafer 
surface are maintained. Consequently, the most effective 

10 CMP systems provide both wafer scale uniformity (i.e., 
uniform material removal over the surface of the 
workpiece) in addition to providing feature scale 
planarity (i.e., removal of small features or surface 
irregularities) . 

15 In a typical prior art CMP system, as shown in 

FIGURE 1, a polishing pad stack includes a relatively 
soft base pad and a relatively stiff upper pad (e.g., a 
polyurethane pad) . The upper pad actually contacts the 
surface of the wafer. The combination of the soft base 

20 pad with the relatively stiff upper pad results in a pad 
stack which is flexible enough to provide uniform 
material removal or wafer scale uniformity across the 
surface of a workpiece, yet stiff enough to smooth out 
the smallest surface irregularities to obtain feature 

25 scale planarity. The slurry used for most CMP systems is 
typically a water based composition which includes . 
suspended colloidal silica particles. 

As may be understood by reference to FIGURES 2A and 
B, CMP is used to planarize the interlevel dielectric 

3 0 material portions of a semiconductor wafer. Typically, 

the surface of the interlevel dielectric materials on an 
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unpolished semiconductor wafer has a pattern on its 
surfaces that results from the dielectric material being 
deposited over a pattern of metal lines. The imbedded 
metal lines within the dielectric material will form the 
electrical connections on the surface of the wafer after 
the dielectric material has been removed by CMP. 

As shown in FIGURES 3A and B, a CMP system may 
alternatively be used to polish away deposited metal 
films such as tungsten or copper by completely removing 
them from a dielectric substrate - except for that 
portion of the deposited layer of tungsten or copper 
which remains in the trenches which are pre-etched into 
the underlying dielectric material . 

The mechanical aspect of the CMP process in the 
removal of small surface irregularities or feature scale 
is the underlying reason why the density of the pattern 
formed by the imbedded metal lines (FIGURE 2) on the 
surface of the wafer plays a big role in the 
effectiveness of a CMP process on removing interlevel 
dielectric materials from the surface of a wafer. In 
other CMP processes, which are designed to remove surface 
metal instead of interlevel dielectric materials, the 
initial planarizat ion of the deposited metal film 
reflects the topography of the underlying dielectric 
material. The removal of metal by a CMP system is more 
dependent on chemistry than on the density of pattern 
lines. The removal of a deposited metal film by CMP has 
a chemical preference that results in a natural "polish 
stop, " wherein the barrier between the resulting main 
conductor and the dielectric material is relatively 
resistant to the slurry on the polishing pads. (The 
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slurry is typically custom blended for the type of metal 
to be removed.) Since polishing is never absolutely 
perfectly uniform across a wafer surface, the deposited 
metal remaining in the trenches etched in the dielectric 
material, from which the conductors are formed, may be 
"dished out" at the top of the trenches while the last of 
the barrier layer of deposited metal film is being 
removed from the surface of the wafer (FIGURE 3) . 
Therefore, for this additional reason, the underlying 
pattern density of the dielectric material plays a 
significant role in the operation of a CMP process to 
remove a deposited metal film. 

The quality of a finished workpiece subjected to a 
CMP process may be expressed in terms of both workpiece 
surface uniformity and workpiece surface planarity. 
Consequently, the particular CMP process parameters 
selected, to include: pad rigidity, slurry composition, 
polish time, relative speed of the wafer and the 
polishing pad, down force imparted on the wafer carrier, 
the dynamics of the workpiece carrier motion, etc., 
assure that the polished workpiece exhibits a desired 
surface uniformity and surface planarity. For example, 
the selection of a CMP process parameter such as the use 
of a comparatively rigid upper pad will tend to planarize 
the surface of a workpiece better than a comparatively 
resilient upper pad. Similarly, the selection of a 
relatively resilient upper pad will provide better global 
uniformity because it can conform to the overall shape 
and contour of the workpiece surface. 

Many prior art systems developed to mathematically 
model a CMP process have focused on calculating the CMP 

AUS01 :223429. 1 



ATTORNEY'S DOWCET PATENT APPLICATION 

090936 . 0395 



system parameters needed to optimize the global 
uniformity of the finished workpiece without considering 
the negative effect that such optimization may have on 
the feature scale planarity of the finished workpiece. 
The results obtained from such prior art CMP modeling 
systems fail to effectively strike an optimized balance 
between both predicting workpiece surface uniformity and 
workpiece surface planarity. 

In the use of CMP processes on interlevel dielectric 
materials, the term "local planarity", which is a feature 
scale measurement, can be used to designate when small 
surface thickness variations on the surface of the 
interlevel dielectric material are within a desired 
tolerance (FIGURE 2 A and 2B) . The analogous measurement 
in the use of a CMP process to remove a deposited metal 
film is the amount of dishing of the deposited metal in 
the top of the trenches in the dielectric material 
(FIGURES 3A and 3B) . In both types of CMP processes, the 
need to achieve local planarity requires that material be 
selectively removed on the feature scale. However, 
selectivity on the feature scale competes, in turn, with 
the need for uniform removal of material on the wafer 
scale. A successful CMP process is one that operates 
effectively in a CMP process regime to ultimately 
optimize both the feature scale and the wafer scale 
measurements on the surface of a polished workpiece. 

It is appropriate to mention that other mathematical 
models directed specifically to feature scale erosion and 
planarization have been developed. One of the first 
published feature scale erosion and planarization models 
appears in the article, "A Two-Dimensional Process Model 
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for Chemical Polish Planarization" by J. Warnock which 
was published in the August 1991 issue of the Journal of 
the Electrochemical Society. An extensive review of Dr. 
Warnock' s model reveals that its parameters were heavily 
5 dependent on the surface pattern of the wafer, thereby 

restricting the applicability of Dr. Warnock ' s model as a 
predictive tool to a narrow range of applications. 

Recently, a closed-form mathematical model 
describing the relationship between pattern density and 

10 the rate of removal of material during the polishing 
process was published by Stine, et al . in an article 
entitled "A Closed-Form Analytic Model for ILD Thickness 
Variation in CMP Processes" which appeared in the 
February 1997 CMP-MIC Conference Journal. The Stine, et 

15 al . model modifies the pressure distribution used in the 
Preston equation for material removal rate based on the 
local pattern density (The Preston equation indicates 
that the material removal rate during a polishing process 
is directly proportional to both the force between the 

20 wafer and the polishing pad, and the relative speed of 
the wafer and the polishing pad.) . 

A feature scale mathematical model based solely on 
pad compression was recently presented by Grillaert, et 
al . entitled "Modeling Step Height Reduction and Local 

25 Removal Rates Based on Pad- Substrate Interactions" at the 
February 1998 CMP-MIC Conference. The Grillaert, et al . 
model represents a superset of the closed- form analytic 
model disclosed in the Stine, et al . article. While 
effective for some applications, the Grillaert, et al . 

30 model still falls short of achieving accurate prediction 
of CMP results because it does not account for the 
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flexural bending of the polishing pads. The flexural 
bending of polishing pads always occurs during an actual 
polishing process . 

The prior art CMP process model disclosed in U.S. 
5 Pat. No. 5,599,423, issued to Parker et al . , is designed 
to simulate and ultimately optimize a semiconductor wafer 
polishing process. The Parker et al . model is based upon 
experimental techniques to optimize the polishing 
parameters. In particular, the Parker et al . model 

10 iteratively varies the process polishing parameters, 
measures the actual polishing results associated with 
each process iteration, and then analyzes the empirical 
data to suggest an optimized polishing process. However, 
the Parker et al . model simply optimizes the global 

15 uniformity of the wafer surface without regard to the 
effect that such global optimization of wafer scale 
uniformity will have on the feature scale planarity of 
the polished wafer. Due to its inherent limitations, the 
Parker et al . model is not capable of optimizing the CMP 

2 0 process parameters in accordance with an initial feature 
scale pattern or in accordance with an intended 
planarizat ion characteristic of the polished wafer. 

Other prior art CMP process modeling or simulation 
techniques lack the capability to calculate a sufficient 

25 number of modeling parameters such that simulation errors 
can be minimized. Such prior art CMP process modeling or 
simulation techniques may produce inaccurate simulation 
results that do not take advantage of the historical 
empirical data associated with a particular set of CMP 

30 process parameters. Furthermore, prior art CMP modeling 
systems are often limited to use with a particular CMP 
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system configuration. Such CMP modeling systems are 
limited by not being capable of processing empirical 
and/or simulated processing results associated with one 
CMP system configuration to model or design another 
5 theoretical CMP system configuration which has a number 
of different physical characteristics and a number of 
different process parameters. 

Finding and expanding a successful mathematical mode 
for a CMP process is a critical need for those involved 

10 in the engineering of manufacturing processes utilizing 

CMP. In the past, expensive experimental techniques have 
been the primary method to define a successful CMP 
process. Therefore, the ability to reliably model and 
thereby predict feature scale surface thickness 

15 variations or wafer scale material removal without having 
to rely on expensive actual experimentation techniques is 
highly desirable. Such a reliable and predictive model 
could be used to define one or more variables in a CMP 
process regime. For example, it is well known that a 

20 high down-force will generally improve the global 

uniformity of the wafer scale of the polished surface by 
providing a more uniform contact between the wafer and 
the polishing pad, but high down forces on the wafer 
carrier also decrease the rate of local planarization of 

25 feature scale because the polishing pad is forced into 

more intimate contact with low areas on the surface being 
polished. Similarly, the use of a comparatively 
resilient polishing pad stack will produce better 
uniformity characteristics but poorer planarization 

30 capability. 
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In addition there is a need for a mathematical model 
for a CMP process which enables the prediction of pad 
life as proper timing of the need for changing polishing 
pads minimizes the cost of a CMP process. 

Accordingly, a robust, predictive model for a CMP 
process is needed that can be used to determine, among 
other CMP process parameters, the optimal wafer carrier 
down force and the best pad rigidity to produce the 
desired surface uniformity and planarity on a workpiece 
for a given CMP process regime. 
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SUMMARY OF THE INVENTION 

The pseudo-physical predictive model for a chemical 
mechanical polishing (CMP) process regime of the present 
invention may be used to predict both the wafer scale 
5 measurement and the feature scale measurement resulting 
after the completion of a CMP process. More 
particularly, the disclosed pseudo-physical CMP model may 
also be used for determining, among other CMP process 
parameters, the optimal wafer carrier down force and the 

10 proper pad rigidity for obtaining the desired amount of 
surface uniformity and planari zat ion . 

The pseudo-physical CMP model disclosed herein is 
based on the premise that the amount of erosion of a 
wafer surface by polishing is proportional to the local 

15 contact force between the polishing pad and the wafer 
being polished. The constant of proportionality is 
termed the "erosion rate coefficient", E, which describes 
the blanket material removal rate as a function of the 
local contact force (e.gr. # microns/second of removal per 

20 pound of force) . The erosion rate coefficient E accounts 
for not only the chemical aspects of a CMP process, but 
also the average amount of wafer to pad contact and the 
relative velocities imposed by the particular CMP process 
regime under study. 

25 The disclosed pseudo-physical mathematical CMP model 

for a CMP process is capable of providing predictions of 
both wafer scale uniformity and feature scale planarity 
on the surface of a finished workpiece. 

Further, the disclosed pseudo-physical CMP process 

3 0 model can generate an optimized set of CMP process 
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parameters that are based on a specified balance between 
wafer scale uniformity and feature scale planarity. 

The present invention also provides a pseudo- 
physical CMP process model that calculates its modeling 
5 parameters in accordance with empirically determined CMP 
results to thereby minimize errors in the CMP simulation 
process . 

The present invention provides a pseudo-physical CMP 
process model that employs interpolation techniques which 

10 enable the effective simulation of CMP process results 
where little or no empirical data exists. 

Finally, the disclosed pseudo-physical CMP process 
model enables the use of simulation data which is 
collected for an existing CMP system. The collected 

15 simulation data aids in the design of a new CMP system 
configuration that has different physical and process 
characteristics than an existing CMP system 
configuration . 
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BRIEF DESCRIPTION OF THE DRAWINGS 

A more complete understanding of the method for 
pseudo-physical modeling of a CMP process, of the present 
invention, may be had by referring to the detailed 
description and claims when considered in connection with 
the drawing FIGURES, wherein: 

FIGURE 1 is a cross-sectional view of the standard 
prior art configuration of a chemical mechanical 
polishing process for the creation of a mathematical 
mode 1 ; 

FIGURE 2A is a cross sectional view of the surface 
of a wafer before CMP having deposited or grown 
interlevel dielectric material which is to be removed by 
CMP; 

FIGURE 2B is a cross sectional view of the surface 
of a wafer after CMP wherein the interlevel dielectric 
has been removed by CMP; 

FIGURE 3A is a cross sectional view of the surface 
of a wafer before CMP having a deposited or grown metal 
film layer which is to be removed by CMP; 

FIGURE 3B is a cross sectional view of the surface 
of a wafer after CMP wherein the deposited metal film 
layer has been removed by CMP; 

FIGURE 4 is a perspective view of a wafer node and a 
pad node about which the mathematical model of the 
present invention is based, which illustrates their 
interaction and the forces therebetween; 

FIGURE 5 is a schematic diagram depicting the spring 
forces used to mathematically model the elastic 
deformation of the pad stack; 
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FIGURE 6 is a schematic diagram of an exemplary 
computer system on which the pseudo-physical CMP process 
model of the present invention may be implemented; 

FIGURE 7 is a flow diagram of the process for 
5 simulating CMP; 

FIGURE 8 is an exemplary wafer scale simulation 
result ; 

FIGURE 8A is an exemplary feature scale simulation 
result ; 

10 FIGURE 9 is a flow diagram of the feature scale 

simulation process ; 

FIGURE 10 is a schematic diagram of a modeled 
dielectric pattern and a corresponding deformation model; 
FIGURE 11 is a flow diagram of the CMP process 
15 parameter optimization process; 

FIGURE 12A is a flow diagram of a simplified CMP 
process model calculation; 

FIGURE 12B is a flow diagram of the CMP process 
model validation; and 
2 0 FIGURE 13 is a flow diagram of the interpolation 

formula generation process. 
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DETAILED DESCRIPTION OF THE INVENTION 
General -Model Applicability 

The pseudo-physical chemical mechanical polishing 
(CMP) process model of a preferred embodiment of the 
5 present invention may be used in conjunction with any 

suitable CMP system configuration designed" for the actual 
removal of material from the surface of a workpiece. 
Although the pseudo-physical CMP process model 10 (FIGURE 
6) and an actual CMP system 12 may be utilized in the 
~5 10 context of any number of different types of workpieces, 

^ the pseudo-physical CMP process model is described herein 

HI in the context of its use for polishing semiconductor 

'f\ wafers. It will also be appreciated by those of ordinary 

W skill in the art that the present invention is not 

15 limited to any particular CMP system configuration or to 
M any specific type of workpieces. Furthermore, because 

Jp the use of a CMP process is generally well known in the 

5 semiconductor fabrication industry, the details of 

semiconductor fabrication will not be described in detail 
20 herein except where necessary for an understanding of the 
present invention . 



Material Re moval Simulation 

Although the present invention may incorporate any 
25 number of suitable modeling techniques, the preferred 
embodiment utilizes methodologies based on the Preston 
equation for material removal during a polishing process: 
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Material Removal Rate {R) = kPS , 
where k (Preston coefficient) is a factor 
representing chemical effects, P is the pressure or force 
imparted by the wafer onto the polishing pad, and S is 
the relative speed between a point on the surface of the 
wafer and a point on the polishing pad. In other words, 
the product PS represents the mechanical effects 
associated with the CMP process. Although the Preston 
equation is relatively simple, the value of the Preston 
10 coefficient k has a significant amount of inherent 

uncertainty because the chemical effects, represented by 
the Preston coefficient k, include the chemical reactions 
between the polishing slurry and the surface of the wafer 
in addition to the availability of the polishing slurry 
15 at the surface of 'the wafer. Consequently, the value of 
the Preston coefficient k can be affected by the chemical 
composition of the polishing slurry, the chemical 
composition of the wafer, the rigidity characteristics of 
the polishing pad, the wafer carrier speed, and the speed 
20 of the polishing table. Indeed, although the Preston 
coefficient k may be held constant to simplify the CMP 
process simulation routine, the value of the Preston 
coefficient k may actually vary across the surface of the 
wafer . 

25 The force P imparted by the wafer on the polishing 

pad may also vary across the surface of the wafer being 
polished. The force P is typically better defined than 
the Preston coefficient k. For example, the average 
value for the pressure or force P can be computed by 

3 0 simply dividing the wafer carrier down force by the 

surface area of the wafer. Additional factors such as 
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the curvature of the wafer backing film and the positions 
of vacuum holes (employed by the wafer carrier to secure 
the wafer) may be analyzed to effectively estimate how 
the force P is distributed across the surface of the 
5 wafer. For the exemplary CMP process application 

described herein, the force P varies quadratically with 
the radius of the wafer. 

The computation of the relative speeds between a 
point on the surface of the wafer and a point on the 

10 polishing pad is accomplished through a straightforward 

application of kinematics. Accordingly, the value of the 
relative speed S between the surface of the wafer and the 
polishing pad can be predicted with near certainty. 
Because the relative speed S does not affect the Preston 

15 coefficient k or the force P, the effect of the relative 
speed S on the CMP process can be predicted with a high 
degree of confidence. 

In the preferred embodiment, the Preston coefficient 
k is one of the pre-determined modeling parameters. 

2 0 Changing the Preston coefficient k results in a 

corresponding scaling of the material removal rate 
associated with the simulated CMP process. In addition 
to the Preston coefficient k, the pseudo-physical 
modeling system 10 of the present invention preferably 
25 employs additional modeling parameters to specify the 

pressure or force distribution across the surface of the 
wafer. As described above, the actual average pressure 
or force across the surface of the wafer can be easily 
calculated if the wafer carrier down force and the 

3 0 surface area of the wafer are known. Thus, the 

pseudo-physical CMP modeling system 10 uses the 
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additional modeling parameters (up to five in the 
preferred embodiment) to approximate the actual pressure 
distribution across the surface of the wafer and to 
simulate the corresponding material removal 
5 characteristics associated with the CMP procedure. It 
should be noted that the pseudo-physical CMP modeling 
system 10 may be configured to hold the Preston 
coefficient k equal to zero when the current 
interrogation or sampling point is located over a groove 

10 or gap in the simulated polishing element or when the 

simulated wafer is overhanging the edge of the simulation 
polishing element. This additional consideration may be 
desirable to compensate for those sampling points where 
the polishing pad does not or cannot contact the surface 

15 of the wafer. 

Computer Requirements 

As shown in FIGURE 6, the effective use of the 
pseudo-physical modeling system 10 is enabled by the use 

20 of a computer 14, supported by a user interface 16, and 
an output device 18. In the preferred embodiment, the 
computer 14 is configured similar to a conventional 
personal computer including suitable processing 
capability 20 and memory capacity 22. 

25 The user interface 16 is configured to enable an 

operator to input data into the pseudo-physical CMP 
modeling system 10 as needed, manipulate the modeling 
results, and otherwise control the operation of the CMP 
process model 10. The output device 18 may be a 

30 conventional computer display terminal, a printer, a 
plotter, or any component suitable for displaying the 
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modeling results and the other information produced by 
the pseudo-physical CMP process model 10. 

Alternatively, the computer 14 may be part of a 
mainframe computing system, part of a computer network 
5 environment, or made an integral part of an actual CMP 
configuration system configuration 12 . Indeed, the 
computer 14 portion of the system may take a variety of 
different forms as long as it includes a sufficient 
amount of processing capability and memory capacity to 

10 support the calculations needed to support the CMP 
process model disclosed below. 

Processor 20 is preferably configured to carry out a 
number of processes (described below) employed in the 
pseudo-physical CMP process model 10. The capability to 

15 carryout such processes may be obtained from software 
programs stored within memory 22 or within a separate 
memory element associated with computer 14 . For the sake 
of convenience, FIGURE 6 depicts a CMP simulation process 
program 24, a model validation process program 28, a CMP 

20 optimization process program 28, an interpolation formula 
generation process program 30, and a feature scale 
simulation process program 31 as discrete functional 
blocks resident within the computer processor 20. The 
processor 2 0 (and/or computer 14) is preferably 

25 configured to interact with a user interface 16, an 

output device 18, and an actual CMP system 12 (via, e.g., 
a CMP controller 32) . 

The computer memory 22 cooperates with the processor 
2 0 in a conventional manner, to store, update, and 

3 0 provide CMP process modeling, CMP process history data, 
and user-defined data to and from the processor 20. 
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Although the computer memory 22 may be arranged in any 
suitable manner, FIGURE 6 depicts a computer memory 22 
organized into a number of distinct databases for the 
sake of clarity. In particular, the pseudo-physical CMP 
process model 10 preferably includes at least one CMP 
process parameter database 34, a modeling parameter 
database 36, an interpolation formula database 38, a 
user-defined CMP data database 40, an empirical CMP 
process history database 42, and a simulated CMP results 
database 44. It should be appreciated that the above 
databases 34, 36, 38, 40, 42, and 44 need not be located 
in a single memory element and that one or more of the 
databases may be stored on a removable storage medium 
such as a floppy disk or a CD-ROM. 

Interface with the Actual CNIP Process 
The actual CMP system configuration in use 12 may 
communicate with the pseudo-physical CMP process model 10 
in a manner that facilitates the real-time optimization 
of the CMP process parameters, termed the "CMP recipe", 
as the surface of a wafer is being polished. To 
accomplish this, the pseudo-physical CMP process model 10 
may provide adjustment instructions to the CMP process 
controller 32 as necessary during the actual running of a 
CMP process. Feedback of empirical measurement data may 
be provided to the pseudo-physical CMP process model 10 
via the CMP controller 32, which may obtain measurement 
data from any number of in- situ workpiece measurement 
systems. Systems for the measurement of semiconductor 
layer thickness, feature scale planarity, and global 
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wafer scale uniformity are well known to those skilled in 
the art and will not be described in detail herein. 

System Operation 
5 The pseudo-physical CMP process model 10 generally 

operates to simulate a CMP process on both a wafer scale 
(for, e.g., uniformity estimates) and a micro- feature 
scale (for, e.g., planarity estimates) . The 
pseudo-physical CMP modeling system 10 is capable of 
7? 10 optimizing the CMP process parameters, e.g., the CMP 

d recipe, in accordance with a user-defined initial surface 

jjji feature pattern on the wafer, a user-defined local 

ft feature profile, and/or a user-defined weighting of the 

Ly relative importance of feature scale planarity to wafer 

7. ■ 15 scale uniformity. 

M The pseudo-physical CMP process model 10 may be 

configured to utilize one or more default initial 
5 modeling parameters for simulating a given CMP process 

simulation (FIGURE 7). The user has the option to alter 

20 such default initial modeling parameters as desired. 

Similarly, the pseudo-physical CMP process model 10 of 
the present invention may provide a number of default 
coefficients for a given modeling parameter to simplify 
the amount of custom data entry required during step 54. 

25 As described in more detail below, step 54 may simply be 
used to input a set of pre-determined initial modeling 
parameters that have been optimized for a particular CMP 
system configuration, a particular wafer size, and/or the 
chemical composition of a particular wafer. 

3 0 Referring now to FIGURE 7, the CMP simulation 

process program 24 employed by the pseudo-physical CMP 
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modeling system 10 is depicted as a flow diagram. The 
CMP simulation process program 24 is carried out by the 
computer 14 . Several user-defined CMP process 
parameters, required to run the CMP simulation process 
5 program 24, are obtained by input to the computer 14 

through user interface 16. In the preferred embodiment, 
an operator is prompted to input the specific data needed 
to run the simulation program 24 . In response to such 
prompts, the user manually inputs the requested data, 
10 selects values generated by the pseudo-physical CMP 
modeling system 10, or accepts one or more default 
values . 

CMIP Process Parameters 
15 As shown in FIGURE 7, the CMP simulation process 24 

preferably includes an initial process parameter step 52, 
during which step 52 the pseudo-physical modeling system 
10 obtains a plurality of CMP process parameters 
associated with a CMP procedure to be performed upon a 
20 workpiece . A "CMP process parameter" is any quantity, 
characteristic, dimension, or other variable that may 
have an effect upon the outcome of an actual CMP process. 
For example, such CMP process parameters may be related 
to the desired CMP recipe, e.g., polish time, speed of 
the polishing element and associated acceleration ramp 
time, workpiece carrier speed and acceleration ramp time, 
workpiece carrier down force acceleration ramp time, 
workpiece carrier sweep range, and workpiece carrier 
sweep speed. A different CMP recipe may also be 
30 specified for distinct polishing stages of the CMP 

process. In addition to polishing recipe information, 
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the CMP process parameters may relate to the physical 
characteristics of the particular CMP system 12 utilized 
during the CMP process, e.g., dimensions and/or rigidity 
of the polishing pad, dimensions of the workpiece to be 
polished, characteristics of the carrier sweep range and 
pivot point, and the like. It should be noted that any 
number of CMP process parameters may be put . into the 
model during step 52 and that the pseudo-physical CMP 
process model 10 may provide one or more default CMP 
process parameters or facilitate the selection of one or 
more preexisting CMP system configurations during step 
52 . The selected CMP process parameters may be stored in 
a database 34 (see FIGURE 7) for future use. 
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Initial Modeling Parameters 

The CMP simulation process program 24 shown in 
FIGURE 7 also involves a second step 54, during which 
step a plurality of initial modeling parameters are put 
5 into the pseudo-physical CMP process model 10. Once 
obtained, the modeling parameters maybe stored in a 
data-base 3 6 and then subsequently accessed by the 
pseudo-physical CMP process model 10. These modeling 
parameters are utilized by the CMP process model 10 to 

10 compute a simulated CMP result associated with, inter 
alia, the CMP process parameters. As explained above, 
each CMP modeling parameter is expressed in terms of a 
quadratic equation. These quadratic equations include a 
number of user-definable coefficients. Thus, step 54 may 

15 obtain or calculate the initial CMP process modeling 

parameters after a number of user-defined coefficients 
have been received. In the context of the present 
invention, each of the initial CMP modeling parameters 
may be a function of one or more of the CMP process 

20 parameters, including: the polishing table speed, the 
workpiece carrier speed, the workpiece carrier down 
force, the polishing element composition or construction, 
the slurry characteristics, and the characteristics of 
the particular film layers being processed, each of which 

25 is associated with the particular CMP procedure to be 
simulated . 

Simulation Control Parameters 
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Referring again to FIGURE 7, the CMP simulation 
process program 24 may also perform a third step 56, 
during which step 56 a number of pre-determined CMP 
process simulation control parameters are put into the 
5 pseudo-physical CMP process model 10. The pre-determined 
simulation control parameters are related to the manner 
in which the CMP process simulation is carried out. For 
example, the pseudo-physical CMP modeling system 10 is 
preferably configured such that a simulated CMP result 

10 includes the data associated with a plurality of discrete 
sampling points on the surface of a given workpiece. The 
use of discrete sampling points on the surface of a given 
workpiece enables the precise comparison of the simulated 
CMP process results to empirical CMP process results 

15 measured at the same points. Accordingly, the 

pre-determined simulation control parameters may be 
related to the actual location and the actual number of 
sampling points on the surface of the wafer, the 
coordinate system used for an output plot, and the time 

20 periods associated with repeated simulations for common 
sampling points. Of course, any number of additional 
settings or variables related to the CMP process 
simulation control parameters may be put in during the 
execution of step 56. 

25 

User Defined Data 

By further reference to FIGURE 7, a fourth step 58 
is preferably performed to input special user-defined CMP 
data associated with the actual CMP process to be 
30 analyzed. The data obtained during step 58 is preferably 
stored in a database 4 0 (see FIGURE 7) for use by the 
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pseudo-physical CMP process model 10. For purposes of 
the present invention, "user-defined CMP data" includes 
theoretical, desired, or actual characteristics of the 
original or the processed wafer that may have an effect 
5 on the execution of the CMP process simulation 24. For 
example, the suggested CMP recipe generated by the 
pseudo-physical CMP process model 10 may be dependent 
upon an initial feature scale pattern on the surface of 
the workpiece, an initial local film thickness profile on 
10 the surface of the workpiece, a desired level of global 
CP wafer scale uniformity on the surface of the workpiece, 

fp or a desired level of feature scale planarizat ion , In 

y ; addition, the CMP simulation procedure conducted by the 

y pseudo-physical CMP process model 10 may be dependent 

15 upon an indicator of the relative importance of wafer 
Q scale uniformity to feature scale planarization for the 

j» surface of the finished workpiece. 

M During the execution of step 58, data indicative of 

an initial pattern on the surface of the wafer and/or an 

20 initial film thickness profile may be put in any suitable 
format. For example, a given pattern may be formed from 
a plurality of arrays, each array being defined by a 
plurality of nodes. Each array may then be defined by 
the number of features contained therein, the length of 

25 its head section, the length of its tail section, the 

number of "hills" and the height of the hills, the length 
of the spaces or "valleys" between the hills, and other 
descriptive elements. In practice, the user may input 
the physical characteristics of the pattern features on 

30 the surface of the workpiece, numerical descriptors, or a 
graphical rendering of the pattern to convey the initial 
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pattern on the surface to be polished to the 
pseudo-physical CMP process model 10. Similarly, the 
initial thickness profile may be rendered in any suitable 
manner, e.g., on a point -by-point basis. 

The performance indicator of the relative importance 
of wafer scale uniformity to feature scale planarity may 
be expressed as a ratio, a percentage, a scaled number, a 
graphical representation, or in any suitable manner. For 
instance, an evenly weighted indicator may cause the 
pseudo-physical CMP process model 10 to optimize both 
wafer scale uniformity and feature scale planarity and to 
generate a set of CMP process parameters or a particular 
CMP recipe to reflect the proper balance between wafer 
scale uniformity and feature scale planarity. However, a 
performance indicator that favors wafer scale uniformity 
may cause the pseudo-physical CMP process model 10 to 
produce an entirely different CMP recipe that is intended 
to increase the global uniformity of the surface of the 
wafer at the expense of feature scale planarization on 
the surface of the wafer. 

Flexibility of Data Input 

It should be noted that inputs associated with all 
of the above steps 52, 54, 56, and 58 of FIGURE 7 need 
25 not be used during the running of the CMP simulation 

process program 24 . Additionally, the inputs associated 
with steps 52, 54, 56, and 58 may be made in a different 
order than that described herein. In addition, the 
relevant data need not be organized, arranged, or 
30 obtained in the specific manner described above. For 
example, the CMP recipe (CMP process parameters) , the 
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initial surface pattern (user defined pattern) , and the 
indicator of global wafer scale uniformity versus feature 
scale planarization (user defined data) may all be lumped 
together as CMP user-defined data. Likewise, the initial 
CMP process modeling parameters and the CMP process 
simulation control parameters may all be lumped together 
as either initial modeling parameters or simulation 
control parameters. In addition, the data and 
information set forth above may be received by the 
computer 14, the memory 22, or the processor 2 0 (or other 
components of computer 14) in any suitable manner and 
such data and information may eventually be stored in the 
memory 22 or routed to the processor 20 for further 
manipulation in accordance with the various processes 
carried out by the computer 14 when executing the 
pseudo-physical CMP process model 10. 

Production of Simulated CMP Result 

After the appropriate input data has been received 
20 by the pseudo-physical CMP modeling system 10, the 
execution of a step 60 in FIGURE 7 causes the 
pseudo-physical CMP process model 10 to perform an 
appropriate modeling routine to obtain a simulated CMP 
process result for the given workpiece . The simulated 
25 CMP process result is eventually stored in a database 44 
(see FIGURE 6) . This data may then be used during 
subsequent runs of the pseudo-physical modeling system 
10. In the preferred embodiment, the simulated CMP 
process result is associated with the CMP process 
3 0 parameters, the initial modeling parameters, and at least 
one element of the user-defined CMP data. In other 
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words, the simulated CMP result generated by the 
pseudo-physical CMP process model 10 is sensitive to a 
change in any of the user-defined CMP data. 

Exemplary Governing Equation 

As described above, the Preston equation relates the 
film thickness at a number of discrete sampling points to 
a number of variables associated with the CMP modeling 
procedure. For example, let (xi y±) be the coordinates of 
N sampling points on the surface of the wafer (where i = 
1 , 2 , 3 , . . . , N) , and let T±, be the film thickness at the 
particular sampling point. Then, the instantaneous film 
thickness for each sampling point on the surface of wafer 
is governed by the following differential equation, which 
is based on the Preston equation: 

^ = -k(t)P(x i ,y i ,t)\\V p (x i ,y i ,t)-V w (x i y i ,t)\\ ; 

where : 

I =1,2,3,.. .tf; 

k(t) = the Preston Coefficient ; 

P - Pressure at the sampling point (xi, yi, t) : 

V p = Velocity of the polishing pad under the sampling 

point (xi, yi, t) ; and 

V w > = Velocity of the wafer at the sampling point 

(xi, yi, t) . 

Erosion Force Relationship 

In practice, the pseudo-physical CMP modeling system 
10 of the present invention performs a number of 
numerical computations to solve the foregoing 
differential equation which is based primarily on the 
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Preston equation. As previously indicated, the relative 
speed of the wafer to the pad, expressed as 

lV p (x i9 y i J)-V w (x n y i J)^ in the foregoing equation, can be 

predicted with near certainty. Therefore, the purpose of 
5 the pseudo-physical CMP process model 10 is to determine 
P, the pressure at the sampling point. Once P is 
determined, the solution to the foregoing differential 
equation is related to the theoretical material removal 
rate associated with the specified CMP procedure. The 

10 theoretical force distribution derived from the various 
modeling parameters used to simulate the CMP process is 
applied during the running of the simulation model to 
predict the film thickness at the specific sampling point 
or points on the surface of the wafer. The change or 

15 decrease in film thickness indicates how much material 

has been removed. It should also be appreciated that the 
rate of material removal from the surface of the wafer or 
the change in film thickness per unit of time, which is 
dependent upon the amount of force applied by the 

20 polishing pad on the surface of the wafer, may differ 

from wafer to wafer and within each wafer. As shown in 
FIGURE 4, the pseudo-physical CMP process model 10 
preferably analyzes the amount of force applied to the 
wafer on a localized nodal scale and determines the 

25 erosion of a particular wafer node in response to the 
localized force at a particular wafer node. 

The nodal contact forces between the polishing pad 
and the wafer to be polished are determined by simulating 
the elastic deformation of the pad stack by first and 

3 0 second abstract mathematical springs, as shown in FIGURE 
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5. The polishing base pad is modeled as a group of 
springs, each having compressibility k ± . The top 
polishing pad is modeled as a group of auxiliary springs 
connected in series to the base pad springs, each having 
5 a compressibility k 2 . The top spring model is 

mathematically tied to the base pad spring model. As 
shown in FIGURE 4, by simultaneously requiring that 
vertical forces be balanced at each point or node on the 
polishing pad and the corresponding point or node on the 

10 wafer to be polished, the shape of the pad, the areas of 
contact, and the local contact force P can be computed. 
Because the properties of the pad and the down- force on 
the surface of the wafer are incorporated directly into 
the model, the effect of changing these physical aspects 

15 of the CMP process can be directly evaluated for any 
particular CMP situation. 

Description of the Pseudo-Physical CMP Process Model 
Algorithm 

2 0 A better understanding of the operation of the 

pseudo-physical CMP process model 10 of the present 
invention for determining the forces P at a pad node 
begins by an understanding of the terminology and naming 
the conventions used. Specifically, W represents a 

25 geometric plane surface that is, on average, parallel to 
the front-side of the wafer (the side to be polished) , 
and C w is a Cartesian coordinate system described with 
respect to plane W by a set of 3 points (x w/ y w , z w ) which 
are established on plane w, such that x w = 0 corresponds 

30 to a point actually on plane W. The wafer surface is 

modeled by dividing it into a collection of points called 
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"nodes" in the three-dimensional space defined in C w . In 
the preferred embodiment, both quad- tree meshing and 
bi-tree meshing techniques are used. 

Each wafer node is numbered, and the individual 
locations of each wafer node are denoted by a subscript. 
For example, the location of the wafer node i is (w w i, y w i, 
z w i) . Associated with each wafer node is a dimension dxi, 
and dyi in the x w and y w directions respectively. The 
dimension dxi, and dyi define a rectangle surrounding the 
node's projection into the z w = 0 plane (i.e., the point 
(w W i, y w i, 0) ) . The defined rectangle on the wafer surface 
is referred to as a "cell" and is given the same number 
as the wafer node it contains. Those cells on the wafer 
surface adjacent to cell i are defined to be those that 
share at least part of one of their cell boundaries with 
that of cell i. The erosion of the features on the wafer 
surface during the polishing process is represented by 
the change in location of the z w coordinate of a nodal 
point on the wafer cell during the CMP process 
simulation. The integer array mi(j), j = 1,2, 3,..., Mi 
contains a list of the Mi cells that are adjacent to that 
of cell i. 

P represents a geometric plane surface that is, on 
average, parallel to the surface of the pad. C p is a 
Cartesian coordinate system whose z-axis is co-linear 
with, but in the opposite direction of C w l s z-axis. The 
coordinates in C p are (x w ±, y w ±, z p ) where when z p = 0 , z p 
corresponds to a point on plane P. The pad surface is 
modeled using a collection of points called "pad nodes" 
in the three-dimensional space defined by C p . In a manner 
identical to the numbering system for the wafer nodes, 
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the pad nodes are numbered and their locations are 
denoted by a subscript. The x-y position of each node on 
the pad surface is the same as the x-y position of the 
counterpart node on the wafer surface. For example, the 
location of a pad node i is (x w ±, y w ±, z p ±) . Associated 
with each pad node is the same spacing in the x w and y w 
directions as used for the wafer node. Thereby each pad 
cell is defined in a manner that is exactly the same as 
the definition for each wafer cell. The deformation of 
the pad's surface caused by the force on the workpiece is 
represented by the change in location of the vertical 
position of the pad nodes during the CMP process 
simulation . 

As shown in FIGURE 4, planes W and P are further 
constrained to be parallel to one another and separated 
from one another by a distance denoted by w 0 such that 
when the vertical position of a node on the pad, z p = w 0 , 
the point on the pad is on the geometric plane IV. Note 
that a location of a wafer node relative to coordinate 
system C p is therefore {x w ±, y w ±, w 0 - z w ±) . 

The following conditions are hereby defined by the 
foregoing system to describe both wafer nodes and pad 
nodes : 

1. ) Wafer Nodes and Pad Nodes Not In Contact: If 

w 0 - z w± > Zpi, then the wafer nodes and the pad 
nodes are not in contact at node i. 

2. ) Wafer Nodes and Pad Nodes In Contact: If w 0 - 

Zwi = z P i, then the wafer nodes and the pad nodes 
are in contact at node i . C denotes the set of 
wafer and pad nodes which are in contact. 
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During a simulation of a CMP process, the x-y 
location of both the pads nodes and the wafer nodes do 
not change. But the vertical locations of both the pad 
nodes z pii and the wafer nodes and z w ±, do change. The 
5 changes in the vertical position z w ± of the wafer nodes 
represent the actual erosion of the wafer during the CMP 
process while the changes in the vertical location z p ± of 
the pad nodes represent the deflection of the pad's 
surface that occurs when the surface of a wafer is 
10 pressed against the pad during the polishing process. 

s - The amount of the deflection of the pad and the erosion 

J3 

yf of the wafer is determined based on the mathematical 

rjj modeling of the forces which are exchanged between each 

pair of wafer and pad nodes at their point of contact . 
§__£_ 15 For each node i in the pad, the pseudo-physical 

nature of the disclosed model assumes that a first 
J5 abstract mathematical spring-like device is connected at 

S the pad node i at the point (x±, y±, 0) . When the 

vertical location of a pad node i z p ± = p 0 , which is 
20 defined as the pad 1 s natural position, the first abstract 
mathematical spring-like device is set to exert no force 
on the pad node i. When the vertical position of a pad 
node z P i is less than p ot the first abstract mathematical 
spring representing the rigidity of the base pad is 
25 referred to as being "compressed." The first abstract 
mathematical spring therefore exerts a force, in the 
positive z p direction, in an amount that is directly 
related to the amount of its compression, p ot -z p ±, and the 
size of the cell, which is designated by the length of 
30 its sides dxi, dyi . When the vertical position of the pad 
node z P i is less than P c , the first abstract mathematical 
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spring representing the rigidity of the base pad is 
referred to as being "stretched" and therefore the spring 
exerts a force in the negative z p direction in an amount 
that is directly related to the amount of its stretching, 
5 Zpi-Zp 0 . The force in the positive z direction exerted on 
the pad node is a first function of the spring deflection 
and the dimensions of the pad node F 2i = F lf (p 0 -Zpit d x , 
dy±) . For clarity, a simple linear model of the abstract 
mathematical spring force of the base pad can be adopted 
10 wherein: 



abstract mathematical spring; 

dxi, dyi is the area of an individual cell. 
For each adjacent pair of nodes, i and j on the top 
pad, there is a second abstract mathematical spring-like 

20 device, representing the rigidity of the top pad, 

connected between the vertical locations z p ± and z P j of the 
adjacent pad nodes i and j. When the vertical position 
of adjacent pad nodes is the same or when z p ± = z P j , which 
is defined as "no bending", the second abstract 

25 mathematical spring exerts no force on either of the two 
adjacent nodes. But when the vertical position of 
adjacent pad nodes is not equal, or z p ± ^ z P j , then the 
second abstract mathematical spring-like device, 
representing the rigidity of the top pad, exerts a force 

30 in the z direction on both adjacent nodes. This force is 
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F 2 i is the spring force at node i; 

rkx is a Hookean spring constant of the base pad; 
Po-Zpi is the amount of compression of the first 
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of a magnitude that is a function of: i) the difference 
between the vertical position of adjacent spring nodes z p ± 
and z P j, ii) the length of l±j of the line segment common 
to cells i and j, and iii) a parameter h, which 
5 represents (but needs not be numerical equal to) the 
thickness of the top polishing pad. The force in the 
positive z direction exerted by the second abstract 
mathematical spring on the pad node i is a second 
function of the difference in the vertical position of 
10 adjacent spring nodes, the length of a common line 

segment, and the thickness of the top polishing pad. 
Therefore it may be denoted by one function 

F 2[i][j] = F 2{ Z P i - 2 pj> l > h )- The spring force F 2[1][J] on the top 
pad attributable to adjacent pad nodes is then defined by 
15 the equation: 



F 2[i][j] ~ rk 2 i z P j z P i Xy* 



wherein rk 2 is a Hookean spring constant; 
2 0 Zpj.z P i is the amount of stretch of the second 

abstract mathematical spring; 

1 is the length of a common line segment 
between cells i and j; 

h is the thickness of the force pad. 

2 5 To determine the total force on a pad node, the two 

previously determined abstract mathematical spring forces 
must be combined. If f± denotes the total force applied 
by an individual pad node i on a corresponding individual 
wafer node i, then when a pad node i and a corresponding 

3 0 wafer node i are not in contact, f± = 0. But when a pad 
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node and a corresponding wafer node are in contact, f± is 
equal to the net force of all of the abstract 
mathematical springs attached to node i. Specifically, 
when the pad node and the corresponding wafer node are 
not in contact, f± = 0 and when in contact, the force at 
the intersection of a pad node and a wafer node is the 
summation of Fii and the total of the F 2 [i][j] forces from 
adjacent nodes as represented below. 

f\ =F \i + S^2[i][jn I .(>)] 

7=1 

Once the total force between a pad node and a wafer 
node has been determined, the erosion of the wafer 
surface or change in film thickness per unit time during 
the CMP process is modeled as the change in the z 
location or the vertical position of the oppositely- 
positioned wafer node. The fraction dz w ±/dt expresses the 
rate at which the height of a node on the surface of the 
wafer erodes. This erosion rate is modeled as a function 
of the net pressure f± on that wafer node by the equation: 

20 dz^ldt^Efi 

where E is a CMP_process erosion rate coefficient 
pre-determined by a comparison to experimental data. 

The overall vertical force of the pad on the wafer 
25 is defined as F p . F p is equal to the sum of all the nodal 
forces fi . Thus, if the parameters of the CMP process 
model are taken to be constant, the overall vertical 
force between the pad and the wafer F p is a function of 
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w 0t the distance between plane W and plane P alone for any- 
given moment in the CMP process simulation. 
Specifically, the further that the wafer is pressed into 
the polishing pad, the higher the force, F p , will be. The 
5 actual CMP process down-force, defined as F , is part of 
the CMP process parameters supplied by the user of the 
algorithm. In typical CMP process applications, a down- 
force on the wafer holding device, F D , is specified by the 
user. Since only a portion of the wafer is represented 
10 in the simulation, the actual process downward pressure, 
F, is taken to be the fraction of the downward force 
borne by the simulated pattern on the surface of the 
wafer and is therefore defined in the algorithm for the 
CMP process to be : 



15 



p _ p pattern 



where A pattern is the area of the simulated pattern and 
A is the total surface area of the wafer. 
20 The CMP process simulation time is broken up into 

several very small time segments, commonly referred to as 
"time steps", each having a duration At. 

The following variables are defined as: 
1 . ) t = time 
25 2 . ) n = number of pad nodes 

3.) tf± n i S h - t\The duration of the polish time 
The overall CMP process simulation algorithm is 
divided into two main parts. They are as follows: 

1.) Erosion Algorithm: The CMP process is taken 
30 through time, the distance between the planes W 
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and P, w 0 is computed at each timestep # and the 
z w coordinate of the wafer node position 
throughout the CMP process is based on the 
nodal forces. 

5 

2.) Contact Force Algorithm: Given a value for the 
distance between planes W and P, w Qf the 
location of each of the pad nodes is first 
computed, the location of each node where the 
10 pad and wafer are in contact is then computed, 

the local force at the contacting nodes is 
computed next, and the sum, F(w 0 ) , of the local 
nodal forces is then computed across the 
surface where the pad contacts the wafer. 
15 The steps used to employ the Erosion Algorithm and 

the Contact Force Algorithm are described below. 

The steps used to employ the Erosion Algorithm are 
as follows: 

1 . Set the start time t = 0 
20 2. Set c, the initial number of pad and wafer 

nodes in contact, c, to empty 

3. Set the end of the first time step to be 

t = t + At 

4. Execute the following steps to determine 
25 the distance between planes W and-P, w Q the force 

between each pad node and each wafer node, f± f and 
the total force of the pad on the wafer, F p (w 0 ) i = 
1 , 2 , 3 , . . . n by : 

30 4.1 Establish two guesses for w Q , the distance 

between planes W and P called w x and w r 
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such that Wi < w Q < w r/ and such that F p 
(w x ) F p (w r ) < 0, i.e., such that the true 
root w Q is "bounded" . 



4.2 Estimate a guess for w a using the 
bisection, secant, Newton, or any other 
number of widely known root finding 
approaches. For example, in the bisection 
method, w Q = (wi + w r j/2. 

4.3 Evaluate F p (w r ) , F p (w x ) , and F p (w Q ) . 

4.4 Estimate a new guess for w± and w r using a 
bisection, secant, or any other number of 
widely known root finding techniques. 

4.5 If the difference between w z and w r is 
sufficiently small or the difference 
between F p (w 0 ) and F is sufficiently small, 
continue to Step 5. Otherwise, repeat 
Step 4.2. 



5. Decrease z w i, i - 1,2, ...,n in the amount 

6. Let the initial number of nodes in contact 
C = C t , the number of nodes in contact at time t. 

7. If t < tfi n i S h resume at Step 3, otherwise 

STOP. 

The steps of the Contact Force Algorithm are as 
follows : 
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1. Set the number of nodes in contact at time 
t, C t = C, the number of initial nodes in contact. 

2 . Establish the following set of equations 

5 to obtain guesses, z p ± for the vertical locations of 

the pad nodes. (Note that forces Fn and F 2 [i][j] are 
purely functions of the pad node locations) : 

The force at the intersection of a pad node and a 
10 wafer node is: 

M x 

F u + E F 2[i][mi0)] = 0, if i e C, 

and the vertical position of pad node i is: 
z P i = w o~ z „i if 1 e c t 

3 . Solve the set of simultaneous algebraic 

15 equations using any number of widely known algebraic 

equation solving techniques. For the example where 
an abstract mathematical model of linear springs is 
used, the Gauss-Seidel (Matrix Solving or Linear 
System solving) method, the Jacobi method, the 

2 0 Conjugate method, the Gradient method,, or the 

Gassian Elimination Method provide suitable 
simultaneous equation solution techniques. 

4. For each node, i = 1, 2, 3, ...,n, if the 
vertical position of a pad node, z p ±> w 0 -z w ±, the 

25 vertical position of a wafer node, then C t —C^Ji t • 

else C, — C5 - {/} 
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5. If the number of nodes in contact at time 
t C t changed during Step 4 , resume at Step 2 . 
Otherwise proceed to Step 6. 

6. Compute the force at the intersection of a 
5 pad node and a wafer node 

fi = 0 otherwise 

7. Compute the total force on the polishing 

10 pad. 

8. STOP. 

During step 60 in FIGURE 7, the CMP process model 10 
preferably obtains at least a wafer scale uniformity 

15 simulation result and a feature scale planarity 

simulation result for the current workpiece . In an 
exemplary embodiment of the present invention, the wafer 
scale uniformity simulation result includes a simulated 
film thickness profile and the feature scale planarity 

2 0 simulation result includes a simulated local pattern 
profile. The simulated film thickness profile may be 
utilized to derive the global wafer scale uniformity 
information, and the local pattern profile may be used to 
derive the local feature scale planarizat ion information. 

25 The global wafer scale uniformity information is related 
to the simulated film thickness calculated at the various 
sampling points, while the feature scale planarity 
information is related to the localized flatness of the 
pattern at particular locations on the wafer surface. 



AUS01 :223429.1 



ATTORNEY'S DOCKET PATENT APPLICATION 

090936 . 0395 

43 



The feature scale simulation step, which may be performed 
by the pseudo-physical CMP process model 10 to obtain the 
feature scale planarity simulation result, is described 
in more detail below. The simulated CMP result may then 
5 contain additional information related to the wafer 

characteristics and the simulated CMP process result may 
be formatted or expressed in any suitable manner for 
maximum utility. 

y 10 System Output 

01 After the modeling routine 60 in FIGURE 7 is 

S completed for the wafer (or after several modeling runs 

EH are completed over a particular time period) , a final 

Ixi step 62 causes the pseudo-physical CMP process model 10 

]\ 15 to produce a suitable output for review by the user. In 

□ the preferred embodiment, step 62 produces an output 

]5 indicative of the wafer scale simulation result and/or 

O; the feature scale simulation result. FIGURE 8A is an 

exemplary wafer scale simulation result (e.g., a film 

2 0 thickness profile) produced by the pseudo-physical CMP 

process model 10. FIGURE 8B is an exemplary feature 
scale simulation result. 

Each plot 70 in FIGURE 8A represents the film 
thickness at certain locations on the surface of the 
25 wafer, e.g., sample points taken along the diameter of 
the wafer. Different plots 70 for the same simulation 
may be associated with the predicted condition of the 
wafer at different processing times. For example, a plot 
72 may represent the film thickness at time t, while 

3 0 another plot 74 may represent the film thickness at a 

future time t+t'. 
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10 



After final step 62 produces one or more outputs 
related to the CMP process modeling results, the CMP 
process simulation program 24 ends. It should be 
appreciated that the CMP process simulation program 24 
may continue with any number of additional tasks and that 
the CMP process simulation program 24 may be incorporated 
into one or more comprehensive processes utilized by the 
pseudo-physical CMP process model 10 or the actual CMP 
system 12 . 



Feature Scale Simulation Output 

FIGURE 8B is an exemplary feature scale simulation 
result (e.g., a local pattern profile) produced by the 
pseudo-physical CMP process model 10. This particular 

15 simulation result includes an initial user-defined 

feature scale pattern 76, e.gr., the original pattern on 
the surface of the wafer put in during step 58 (see 
FIGURE 7) or a derivative thereof. The CMP process 
simulation result may also include one or more plots 78 

20 representing various stages during the simulated CMP 
process. Such CMP process simulation results may be 
utilized to determine the effect that CMP process 
parameters have on the local planarization of the surface 
of the wafer. 

25 FIGURE 9 is a flow chart which depicts the feature 

scale simulation process 31. The feature scale 
simulation process 31 may be performed by the 
pseudo-physical CMP process model 10 during the execution 
of the CMP simulation process program 24 . The feature 

30 scale simulation process 31 begins with a step 150, which 
causes the pseudo-physical CMP process model 10 to obtain 
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the initial feature scale pattern associated with the 
surface of the workpiece . As described above in 
connection with step 58 (see FIGURE 7) , the initial 
feature scale pattern may be stored in the computer 
5 memory 22 for subsequent access during the execution of 
the feature scale simulation process 31. Referring now 
to FIGURE 10, an exemplary feature scale pattern 170 is 
illustrated. This feature scale pattern 170 is shown as 
if the wafer is face-down in a position ready for 

10 polishing by a polishing pad stack (see FIGURE 1) . As 
shown, the feature scale pattern 170 is preferably 
represented by a plurality of nodes 172 which are 
connected by line segments for the purpose of simulating 
the CMP process in accordance with the present invention. 

15 The spacing between the nodes 172 may be selected in any 
suitable manner, e.g., to provide a sufficiently accurate 
model of the feature scale pattern 170. 

Pad Deformation Model Output 

2 0 Referring back to FIGURE 9, a step 152 is performed 

to apply the mathematical model of the deformation of the 
polishing element, e.g., the polishing pad stack used by 
the actual CMP system 12 . The preferred model of pad 
deformation is based upon the physical characteristics of 

25 the polishing pad acting like abstract mathematical 
springs connected together in series (FIGURE 5) . An 
exemplary mathematical model 174 of a deformation pad is 
shown in FIGURE 10 where the modeled pad, when in its 
deformed state, closely simulates the characteristics of 

30 a polishing element used in an actual CMP system 12. In 
the preferred embodiment, the mathematical model 174 of 
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the polishing element is defined at individual pad nodes 
176 that correspond (relative to the horizontal axis) to 
individual wafer nodes 172. Further, the mathematical 
pad deformation model 174 assumes that adjacent nodes 176 
on the pad are connected together by line segments. 

At any given time during the CMP process simulation, 
the feature scale pattern 170 on the surface of the wafer 
is considered to be rigid (not deformed) and the 
polishing pad is considered to be deformable in 
accordance with the mathematical model 174. The 
mathematical model of pack deformation 174 is preferably 
defined by a plurality of primary force/displacement 
elements 178 (abstract mathematical springs) that are 
capable of compressing in response to a load. In the 
exemplary embodiment shown in FIGURE 10, each primary 
force/displacement element 178 in the pad is associated 
with one of the pad nodes 176. Each primary 
force/displacement element 178 in the pad may then be 
"linked" or otherwise associated with at least one 
additional adjacent primary force/displacement pad 
element 178. In the preferred embodiment, the primary 
force/displacement elements 178 in the pad are then 
"coupled" to each of their adjacent primary 
force/displacement elements 178 via secondary 
force/displacement elements 180. While any number of 
modeling parameters may be employed to characterize the 
primary and secondary force/displacement pad elements 
178, 180 -- in the preferred embodiment the displacement 
pad elements 178, 180 are physically modeled to function 
as linear springs. 
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After step 152 (FIGURE 9) acquires the pad 
deformation model 174, a step 154 is preferably performed 
to initialize a displacement 182 of the wafer relative to 
the polishing pad. In practice, the displacement 182 of 
5 the wafer may relate to an amount of downward travel of 
the wafer carrier during the CMP procedure. The initial 
static displacement of the wafer with respect to the 
polishing pad stack may be set at any suitable value. In 
response to this initial static wafer displacement, a 
10 step 156 causes the pseudo-physical CMP process model 10 

CP to determine (by simulation) the deformation of the 

=£1 

i a polishing element. During the execution of step 156, the 

pseudo-physical CMP process model 10 may suitably 
yj generate a simulated contact profile between the wafer 

!\ 15 and the pad and/or a localized force profile associated 

Q with the pad deformation model 174 in relation to the 

"« current state of the feature scale pattern 170 on the 

M wafer. The contact profile between the surface of the 

wafer and the pad may identify which nodes 172 of the 
20 feature scale pattern 170 on the surface of the wafer are 
in contact with the corresponding nodes 176 of the pad 
deformation model 174. As shown in FIGURE 10, some parts 
of the simulated pad may not be in contact with 
corresponding parts of the wafer, and the secondary 
25 force/deflection elements 180 may limit the extension of 
some of the primary force/deflection elements 178 that 
would otherwise be extended to contact the simulated 
pattern on the surface of the wafer. The force profile 
may then identify the amount of localized force that is 
30 imposed upon the nodes 172 of the feature scale pattern 
on the surface of the wafer by the nodes 176 of the pad 
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deformation model by analyzing the appropriate pad 
force/deflection elements 178, 180 with respect to the 
pad displacement 182. 

Contact Profile Computation 

The contact profile between the wafer and the 
polishing pad is preferably computed by utilizing a 
static force equilibrium equation at each node of the pad 
deformation model 176. For any given pad node, the 
downward displacement of its primary force/deflection 
element induces a corresponding upward force, just like a 
linear spring. The greater the displacement -- the 
greater the force. The relative displacement of this pad 
node with respect to its adjacent pad nodes induces 
additional forces (either upward or downward, depending 
on the relative position of the pad node to its adjacent 
node) . The pseudo-physical CMP process model 10 
preferably assumes that the summation of the nodal forces 
for any given nodal contact between the pad and the wafer 
is equal to zero. When summation of the nodal forces 
equaling zero is imposed mathematically for any given pad 
node, an equation in terms of the displacement of the 
given pad node and the displacement of those pad nodes 
adjacent to the given pad node is formed. Combining such 
equations for each one of the pad nodes 176 of the pad 
deformation model produces a set of simultaneous 
algebraic equations which mathematically represents the 
displacement of the pad nodes 176 in the pad deformation 
model. In that set of simultaneous algebraic equations, 
the force-balance equation for any node on the polishing 
pad that has previously been determined to be in contact 
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with a corresponding node on the wafer is replaced with 
an equation holding that the vertical position of a node 
on the polishing pad is equal to the vertical position of 
a corresponding node on the wafer at that point. 
5 After the set of simultaneous algebraic equations 

has been solved, the displacement of each node on the pad 
is examined. If the position of a node on the polishing 
pad is determined to be higher than that of the position 
of a corresponding node in the wafer at that point, then 

10 the node position in the polishing pad is automatically 
set to be equal to the node position of the wafer and 
this node is identified as a "possible contact node." 
The set of simultaneous algebraic equations is then 
solved once again -- this time replacing the force 

15 balance equation for each "possible contact node" with a 
specified pad displacement relative to plane P. The 
specified pad displacement is the position of the wafer 
surface relative to plane W at that node. Similarly, if 
the vertical position of a node on the polishing pad, 

2 0 previously thought to be a "possible contact node, " and 

then later is determined not to be in contact with a 
corresponding node on the wafer, then that pad node is 
removed from the list of "possible contact nodes" and its 
corresponding position equation is replaced with the 
25 force balance equation. Solving the set of simultaneous 
algebraic equations and then subsequently determining 
"possible contact nodes" is repeated until no changes 
occur between iterations. At that point, the "possible 
contact nodes" of the polishing pad are determined to be 

3 0 in contact with the corresponding wafer nodes and are 

identified as "contact nodes." The force imposed on the 
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wafer nodes by each contact node of the polishing pad is 
then computed using the force relationships associated 
with the primary and the secondary force/deflection 
elements 178, 180 at each pad node. 

5 

Updated Displacement Model 

After step 156 generates an appropriate force 
profile for the current displacement 182, a quaere step 
158 is then performed. The quaere step 158 compares the 

10 sum of the node forces contained in the localized force 
profile to the current (or simulated) down force 
associated with the particular area of the wafer being 
analyzed. This current (or simulated) down force may be 
a user-defined quantity or it may be derived from a known 

15 carrier down force associated with the entire wafer and 
the known area of the wafer under analysis. If the 
quaere step 158 determines that the sum of the nodal 
forces does not substantially equal the local down force 
(or, alternatively, if the difference between the sum of 

2 0 the nodal forces and the current (or simulated) down 

force does not fall within a desired tolerance) , then a 
subsequent step 160 is performed to adjust the pad 
displacement 182 by a suitable amount. The adjustment of 
the pad displacement 160 causes a change in the pad 

2 5 deformation model 174, the current contact profile, and 
the current localized force profile. Accordingly, step 
156 and the quaere step 158 are preferably repeated for 
the updated pad displacement 182 . 

Steps 156, 158, and 160 preferably form a processing 

30 loop that causes the pad displacement 182 to be adjusted 
in a suitable manner until the quaere step 158 determines 
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that the sum of the node forces on the polishing pad is 
approximately equal to the local down force applied to 
the wafer. If these forces substantially balance, then a 
step 162 is performed to cause the pseudo-physical CMP 
5 process model 10 to simulate the erosion of the wafer for 
a given time period using the differential equation for 
dTi/di shown above. In the preferred embodiment, the 
erosion of the wafer surface is simulated by adjusting 
the positions of the appropriate nodes 172 of the feature 

10 scale pattern on the wafer surface by an amount that may 
be dependent upon the current CMP process parameters, 
such as: the slurry composition, the wafer composition, 
and the like. Thus, for a given time, soon after CMP 
processing begins, a local pattern on the wafer surface 

15 may closely resemble the original wafer surface pattern 
(see FIGURE 10) . At a later time, the pseudo-physical 
CMP modeling system 10 may leverage historical simulation 
data to enable the simulation of the erosion of the wafer 
surface during the CMP process over time. For example, a 

20 subsequent iteration of the feature scale simulation 

process 31 may utilize a partially planarized version of 
the wafer surface pattern rather than the initial wafer 
surface pattern described above in connection with step 
150 (FIGURE 9) . 

25 

Optimization Process 

FIGURE 11 depicts the CMP optimization process 
program 2 8 that may be performed by the pseudo -physical 
CMP process model 10 to generate a preferred set of CMP 
3 0 process parameters for use during an actual CMP 

procedure. The CMP process parameters are preferably 
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optimized to produce a CMP recipe in response to a 
desired CMP result. The optimization process 28 may be 
carried out in conjunction with the CMP process 24 or as 
a separate and distinct process. 
5 The CMP optimization process program 28 preferably 

begins with an initial step 86 which causes the 
pseudo-physical CMP process model 10 to retrieve a 
current simulated CMP result from, e.g., database 44 
located in the memory 22 (see FIGURE 6) . Following step 

10 86, another step 88 is performed to retrieve the 

optimization parameters for the CMP optimization process 
28. These optimization parameters may be retrieved from 
the database 40. In the context of this description, an 
optimization parameter may be any characteristic, 

15 quantity, or feature associated with an "ideal" wafer as 
processed by the actual CMP system 12 . The preferred 
embodiment of the pseudo-physical CMP process model 10 
may include any number of optimization parameters related 
to: the initial feature scale pattern; the initial film 

20 thickness profile; the relative importance of wafer 

uniformity versus wafer planarization; or the intended 
planarization or uniformity result. All of these 
parameters may be user-defined. 

Following the completion of a step 88, a subsequent 

25 step 90 is performed to cause the pseudo-physical CMP 
process model 10 to compare the characteristics of the 
current simulated CMP result to the corresponding 
characteristics indicated by the optimization parameters. 
For example, in the context of feature scale planarity, 

3 0 step 90 may analyze an error between the intended pattern 
on the surface of the wafer and the simulated pattern, 
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which error is contained in the current simulated CMP 
results. Similarly, with respect to wafer scale 
uniformity, step 90 may analyze the difference between 
the intended film thickness and the simulated film 
thickness measurements . 

A quaere step 92 is preferably performed to 
determine whether the error between the simulated CMP 
result and the intended CMP result has been substantially 
minimized. The CMP optimization process 28 may utilize 
any number of techniques to analyze the simulation error, 
e.gr., curve fitting, least-squares, averaging, and the 
like. If the quaere step 92 determines that the 
simulation error is substantially minimized, then the 
current simulation is considered acceptable and the 
quaere step 92 may prompt another step 94, which both 
saves and displays the current CMP process parameters. 
Following the completion of step 94, the CMP optimization 
process 28 ends. The operator may then apply the 
optimized CMP process parameters to the actual CMP 
procedure. If the quaere step 92 finds that the 
difference between the simulated and the intended CMP 
results have not been minimized, then a step 96 may be 
performed to cause the pseudo-physical CMP modeling 
system 10 to adjust at least one of the CMP process 
parameters. For example, step 96 may vary the polish 
table speed, the wafer carrier down force, the polish 
time, the dimensions or sweep range of the carrier, or 
the like. Step 96 may also hold one or more CMP process 
parameters fixed (in response to a user input or 
automatically) to simplify subsequent processing. 
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After step 96 adjusts the initial CMP process 
parameters, step 98 causes the pseudo-physical CMP 
modeling system 10 to perform the modeling routine with 
the updated CMP process parameters in place. The 
5 modeling routine was described above in connection with 
FIGURE 9. Steps 90 and 92 are repeated to analyze the 
new simulated CMP result. In this manner, steps 90, 92, 
96, and 98 form a processing loop during which one or 
more of the CMP process parameters are optimized in 

10 accordance with the intended CMP results. The particular 
set of CMP process parameters may be saved in a database 
34 (see FIGURE 6) for future reference or to initialize a 
subsequent CMP process optimization routine. 

Those skilled in the art will appreciate that the 

15 CMP optimization process 28 (or a modified version 

thereof) may be employed to optimize those CMP process 
parameters related to a theoretical CMP system. Such 
simulations may facilitate the design and development of 
new CMP systems without the cost and labor associated 

2 0 with actual experimentation and prototyping. For 
example, after a particular CMP procedure has been 
optimized for an existing CMP system, the CMP 
optimization process 2 8 may be performed to vary at least 
one of the optimized CMP process parameters to thereby 

25 define an updated CMP process parameter set. Then, a 
second simulated CMP result may be obtained using the 
updated CMP process parameters. In this manner, an 
operator can efficiently simulate and optimize the 
performance of a new CMP system prior to building an 

30 actual prototype. 
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It should be noted that the CMP optimization process 
28 is not limited to the particular optimization protocol 
described above. Indeed, the CMP optimization process 28 
may utilize any suitable optimization technique to 
accomplish the same results. For example, the CMP 
optimization process 28 may utilize minimization or 
maximization methodologies to optimize one or more values 
associated with quantifiable CMP results. In the 
preferred embodiment, the CMP optimization process 28 may 
endeavor to maximize the quality of the CMP result by 
adjusting the CMP process parameters in accordance with 
the following technique. A quality measurement (Q) is 
defined as follows: (Q) = ( ) (uniformity) + (1- ) 
(planarity) , where the user selects a value for (the 
relative importance of uniformity compared to planarity) 
between zero and one. Thus, the CMP optimization process 
28 may iteratively adjust the CMP process parameters to 
maximize Q, subject to the particular value of 

2 0 Model Validation Process 

In the preferred embodiment, the pseudo-physical CMP 
modeling system 10 is capable of validating one or more 
of the processing parameters associated with a simulated 
CMP result such that the error between the simulated CMP 

25 result and the empirical CMP result (using the same CMP 
process parameters) is substantially minimized. Such 
model validation is desirable to determine whether a 
given CMP model can reproduce experimental data through 
manipulation of one or more of its modeling parameters, 

30 to establish rules for computing the modeling parameters 
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for new configurations, and to determine the reliability 

of such rules and the accuracy of the CMP simulations. 

FIGURE 12A is a simplified flow diagram illustrating 

that a partially polished wafer is measured to validate 
5 the mathematical model and obtain updated CMP process 

settings to be plugged back into the CMP process model to 

refine its output. 

FIGURE 12B is a detailed flow diagram of the model 

validation process 2 6 that may be performed by the 
T3 10 pseudo-physical CMP process model 10. Generally, the 

5} model validation process 2 6 enables the pseudo-physical 

yff CMP process model 10 to validate its CMP simulation 

jf: results by comparing the simulated results to 

UJ corresponding empirical results and then adjusting one or 

1^ 15 more of the modeling parameters to increase the accuracy 

M of the simulation. Thus, the pseudo-physical CMP process 

Jc model 10 may be configured to leverage historical 

empirical data to improve the overall performance of the 

CMP process simulations. The model validation 26 may be 

2 0 performed for any number of modeling parameters and for 

different CMP process parameters to enable accurate 
simulations for a wide variety of different CMP 
processing recipes . 

The model validation process 26 begins with step 
25 112, which causes the pseudo-physical CMP process model 
10 to retrieve a particular or current CMP simulation. 
Step 112 is similar to step 86 described above in 
connection with FIGURE 11. The CMP simulation retrieved 
during the execution of step 112 is generated by the 

3 0 pseudo-physical CMP process model 10 in accordance with 

the present invention. An actual run of a CMP procedure 
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is conducted during the execution of step 114. This CMP 
procedure is performed by the actual CMP system 12 in 
accordance with the CMP process parameters designated 
during the corresponding CMP simulation. Step 114 may be 
performed in response to control signals from the pseudo- 
physical CMP process model 10 or in response to operator 
inputs to the actual CMP system 12. In addition, step 
114 may be performed in an interactive manner as the 
pseudo-physical CMP process model 10 produces a 
substantially real-time CMP simulation. 

Following step 114 a step 116 may be performed to 
measure an empirical CMP process result associated with 
the wafer polished during the execution of the actual CMP 
procedure. The empirical CMP process result is 
preferably stored in a database 42 (see FIGURE 6) for 
subsequent use by the pseudo-physical process model 10. 
In the preferred embodiment, step 116 may be performed by 
a measurement system incorporated into the actual CMP 
system 12 and the pseudo-physical CMP process model 10, 
or by one or more separate measurement systems . The 
actual CMP system 12 may include any number of in- situ 
wafer surface measurement devices to facilitate 
substantially real-time optimization of one or more of 
the CMP process parameters in response to the current 
simulation results. Such measurement devices (and other 
suitable wafer surface measurement systems) may be 
employed for purposes of executing step 116. Such 
surface measurement systems and devices are well known to 
those skilled in the art and accordingly are not 
described in detail herein. 
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Step 116 preferably measures at least the global 
wafer scale uniformity of the processed wafer (derived 
from a film thickness profile) and the local feature 
scale planarity of the processed wafer (derived from a 
5 local feature pattern profile) . Thus, step 116 obtains 
an empirical CMP result that contains a wafer scale 
empirical CMP result and a feature scale empirical CMP 
result. The global wafer scale uniformity may be derived 
by measuring the film thickness at a number of points on 

10 the surface of the wafer. In an exemplary embodiment, 
the film thickness profile is measured with a thin-film 
thickness measurement system, such as Optiprobe system. 
The local planarity may be derived by scanning the 
surface of the wafer and analyzing the small scale 

15 features. For example, the preferred embodiment utilizes 
a surface profile measurement system which is 
commercially available as a stand alone system. Other 
suitable techniques may be utilized during step 116 to 
measure uniformity, planarity, or other characteristics 

20 of the wafer, e.g., systems that employ reflective or 
refractive optics or systems that employ micrometer or 
observational techniques. For the sake of compatibility, 
the measurement points preferably correspond to the 
sampling points used by the pseudo-physical CMP process 

25 model 10 to produce the current simulated CMP result. 

After step 116 obtains the requisite amount of 
empirical CMP data, a step 118 is performed to compare 
the simulated CMP result to the empirical CMP result. In 
practice, the pseudo-physical CMP process model 10 may 

30 compare the simulated and empirical film thickness 

profiles and the simulated and empirical local feature 
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profiles. During step 118, the actual wafer scale 
empirical CMP result and the actual feature scale 
empirical CMP result are compared to their respective 
simulated counterparts. The pseudo-physical CMP process 
5 model 10 may use any suitable comparison technique during 
task 118 to compare the empirical and simulated results. 
For example, the various sampling points may be analyzed 
on an individual basis or a plurality of sampling points 
associated with a given measurement may be processed in a 

10 collective manner. Alternatively, step 118 may employ 
any number of conventional curve fitting techniques, 
least-squares techniques, or the like. 

In accordance with an exemplary embodiment of the 
CMP process model 10, a quaere step 120 is performed 

15 after the empirical and simulated CMP results have been 
obtained. The quaere step 12 0 determines whether a 
simulation error (which may be obtained during step 118) 
is substantially minimized. Such a simulation error may 
be determined for individual sampling points or for a 

2 0 collected or averaged quantity with a number of sampling 
points. For example, in the preferred embodiment, the 
difference is mathematically determined by summing the 
individual differences between the simulated and 
empirical results at various data points. The sum of 

2 5 these individual differences is dependent upon the values 

of the modeling parameters used by the pseudo-physical 
CMP process model 10 . 

If the quaere step 12 0 determines that the current 
simulation error has been minimized, then a step 122 is 

3 0 performed to store the current CMP process parameters for 

use with a subsequent CMP modeling routine. In other 
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words, the optimized CMP process parameter set may be 
used to conduct future simulations in a confident manner. 
The optimized CMP process parameters may be stored in a 
database 36 of the memory 22 (see FIGURE 6) . After the 
completion of step 122, the model validation process 26 
ends. It will be appreciated by those of ordinary skill 
in the art that the model validation process 26 may be 
performed in conjunction with the CMP optimization 
process 28 (FIGURE 11) in an iterative or combined manner 
to obtain optimized CMP process parameters for a given 
CMP procedure . 

If quaere step 120 determines that the current 
simulation error is unacceptable, then the model 
validation process 26 leads to a step 124. Step 124 
causes the pseudo-physical CMP process model 10 to adjust 
at least one CMP process parameter to obtain an updated 
CMP process parameter set . The pseudo-physical CMP 
process model 10 may be suitably configured to 
systematically adjust the CMP process parameters in a 
number of ways and in any order. Step 124 may also cause 
the pseudo-physical CMP process model 10 to hold at least 
one of the CMP process parameters fixed during the 
adjustment procedure to facilitate efficient and speedy 
optimization. The specific CMP process parameters to be 
fixed may be designated by the user prior to the model 
validation process 26, during the model validation 
process 26, or by default. It should be noted that, 
because the simulation error is dependent upon the 
various CMP process parameters, the model validation 
process 26 is reduced to a multi-variant optimization 
problem. 
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Following step 124, a step 126 preferably causes the 
pseudo-physical CMP process model 10 to perform a 
subsequent modeling routine to obtain an updated 
simulated CMP result associated with the updated CMP 
process parameters. Step 126 is similar to step 98, 
described above, in connection with FIGURE 11. After 
step 12 6 obtains the new simulated CMP result, the model 
validation process 26 is re-entered at step 118 to re- 
compare the empirical CMP result with the updated 
simulated CMP result. Thus, steps 118, 120, 124, and 126 
are preferably repeated until the error between the 
simulated CMP result and the empirical CMP result has 
been substantially minimized. In other words, the 
simulated CMP result is altered until a "best fit" 
relative to the empirical CMP result has been obtained. 
This processing loop causes the pseudo-physical CMP 
process model 10 to self -optimize its CMP process 
parameters such that the simulated CMP result 
substantially matches the empirical CMP result. Although 
the present invention employs non- linear regression 
techniques to optimize the CMP process parameters, any 
number of suitable methodologies may be utilized to 
determine the particular CMP process parameters used for 
a given simulation. 

Interpolation Formula Generation Process 
FIGURE 13 is a flow diagram of the interpolation 
formula generation process 3 0 performed by an exemplary 
embodiment of the pseudo-physical CMP process model 10. 
The interpolation formula generation process 3 0 
preferably begins with a step 132, during which step the 
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pseudo-physical CMP process model 10 obtains a plurality 
of optimized CMP process parameter sets, each set being 
associated with a pre-determined set of CMP process 
parameters. In other words, each of the optimized CMP 
5 process parameter sets is configured for use with the CMP 
modeling routine to obtain a simulated CMP result for a 
wafer processed during the CMP procedure using a set of 
pre-determined CMP process parameters. The optimized CMP 
process parameters are preferably obtained during the 

10 model validation process 26 (see FIGURE 12B) or by an 
equivalent CMP process parameter estimation technique. 
Task 132 may obtain the optimized CMP process parameters 
from database 3 6 (see FIGURE 6) . 

Step 132 leads to a step 134, which causes the 

15 pseudo-physical CMP process model 10 to develop a 

plurality of interpolation formulas associated with the 
optimized CMP process parameter sets. The interpolation 
formulas may be utilized by the pseudo-physical CMP 
process model 10 to predict CMP process parameter values 

20 for new and untested sets of various CMP process 

parameters. Those skilled in the art will appreciate 
that any number of methods may be used to form the 
interpolation rules for the CMP process parameters. The 
most straightforward manner is to simply express each CMP 

25 process parameter as a linear or as a quadratic function. 
Such a methodology may be acceptable when a large amount 
of empirical data is accumulated to encompass an adequate 
number of CMP procedures. However, this procedure may 
not be adequate when a limited amount of empirical data 

30 exists. 
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The preferred embodiment endeavors to generate 
interpolation formulas associated with those CMP process 
parameters that dictate performance and/or provide for 
reliable interpolation rules for the other CMP process 
parameters. One technique that accomplishes this goal 
limits the number of variable CMP process parameters to 
the following: 

T = table speed 

C = carrier speed 

F = carrier down force 
The pseudo-physical CMP process model 10 uses a 
combination of these parameters for interpolating the CMP 
process parameters; this combination is defined as II, 
where : 

II = r e » C e2 (T-C) e3 F e * 

If x 2 is an objective function of a minimization 
problem and x 2 = (ei, e 2 , e 3 , e 4 ) then an interpolation 
rule may be obtained by calculating the different 
possible combinations and analyzing x 2 . As an 
approximation to this approach, the exponents can be 
restricted to the values -1, 0, or 1 (with e 3 ^ 1) . The 
best overall "fit" is then selected for purposes of this 
interpolation formula . 

After the interpolation formulas are adequately 
developed in step 134, a step 136 may be performed such 
that the pseudo-physical CMP process model 10 receives 
CMP data associated with a particular CMP procedure. As 
described above, this CMP data may be related to the CMP 
process parameters and/or the intended CMP results (e.g., 
uniformity versus planarity balance, wafer pattern, or 
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the like) . Step 136 may be performed in response to 
user-defined inputs or in accordance with other 
instructions received by the pseudo-physical CMP process 
model 10. 

5 A step 138 is preferably performed in response to 

the CMP data received in step 136. Step 138 causes the 
pseudo-physical CMP process model 10 to generate an 
interpolated modeling parameter based upon the current 
CMP data. Step 138 utilizes one or more interpolation 
2 10 formulas to generate the interpolated modeling parameters 

4; for the given CMP data. Following step 13 8, the 

yj interpolated formula generation process 3 0 may end. The 

Jf! interpolated CMP process parameters may then be stored or 

Uj immediately used to obtain a simulated CMP result for the 

l^, 15 specified CMP procedure. 

^ As depicted by the ellipses and the connecting arrow 

J: in FIGURE 13, a step 140 (which may be performed at a 

2 later time or in connection with a separate process) may 

be performed to produce a plurality of CMP process 
20 parameters in response to an intended CMP result, where 
the CMP process parameters are obtained through the 
optimization procedures described above, which procedures 
involve the use of the interpolation formulas. Thus, the 
interpolated formula generation process 30 (or a related 
25 process) maybe performed by the pseudo-physical CMP 
process model 10 to obtain suggested CMP process 
parameters for an untested CMP procedure if an intended 
CMP result is known beforehand. 



30 
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Utilization of CMP Process Parameters 
The CMP process parameters obtained during step 14 0 
may be subsequently (or substantially concurrently) 
applied to the actual CMP system 12, as depicted in a 
step 142 related to the interpolated formula generation 
process 30. Step 142 may cause the pseudo-physical CMP 
process model 10 to communicate the CMP process 
parameters to the CMP controller 32 with the actual CMP 
system 12 (see FIGURE 6) . Alternatively, an operator may 
make a record of the CMP process parameters and adjust 
the actual CMP system 12 in an appropriate manner to 
effect the desired settings for CMP process parameters. 
Following step 142, a step 144 causes the actual CMP 
system 12 to perform a CMP procedure in accordance with 
the set of CMP process parameters produced during step 
140. Step 144 is similar to step 114 described above in 
connection with FIGURE 12B. Following step 144, the 
process 3 0 ends. As described above, the empirical CMP 
results obtained during step 144 may be used during the 
model validation process 26 (FIGURE 12B) . 

Other Embodiments 

In summary, the present invention provides an 
improved CMP mathematical model for a CMP process model 
25 that is capable of providing both wafer scale uniformity 
simulations and feature scale planarity simulations. The 
disclosed CMP modeling system can generate an optimized 
set of CMP process parameters based on a specified 
balance between wafer uniformity and planarity. In 
30 addition, the pseudo-physical CMP process model computes 
its modeling coefficients in accordance with empirically 
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determined CMP results to thereby reduce simulation 
errors. The pseudo-physical CMP process model employs 
interpolation techniques to its modeling parameters to 
effectively simulate CMP results for which little or no 
empirical data exists, and it can process simulation data 
collected for an existing CMP system to assist in the 
design of a new CMP system having different physical and 
processing characteristics than the existing CMP system. 

The present invention has been described above with 
reference to preferred exemplary embodiments. However, 
those skilled in the art will recognize that changes and 
modifications may be made to the preferred embodiment 
without departing from the scope of the present 
invention. For example, the present invention is not 
limited to the particular modeling, optimization, or 
interpolation techniques described herein. In addition, 
the various hardware components may differ than that 
shown and described herein and the various processes need 
not be performed in the precise manner described therein. 
These and other changes or modifications are intended to 
be included within the scope of the present invention, as 
expressed in the following claims. 
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